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In supersymmetric generalizations of the Standard Model, all stable Q-balls are associated with 
some flat directions. We show that, if the flat direction has both the baryon number and the lepton 
number, the scalar field inside the Q-ball can deviate slightly from the flat direction in the ground 
state. We identify the true ground states of such nontopological solitons, including the electrically 
neutral and electrically charged Q-balls. 
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' Nontopological solitons, Q-balls appear in theories with some scalar fields carrying a global conserved quantum 

1 . J number when the scalar potential meets the conditions for Q-ball stability with respect to a decay into scalar parti- 
O • cles [1-3]. These conditions are satisfied in supersymmetric generalizations of the Standard Model for the scalar fields 

' carrying the baryon and lepton numbers [4] . In theories with gauge-mediated supersymmetry breaking, Q-balls can be 
^ ■ stable with respect to decay into both scalar particles and fermions if the vacuum expectation value (VEV) inside the 
Q-ball corresponds to a flat direction [5]. These stable objects can form in the early universe from the fragmentation 
I— i, of the Affleck-Dine [6-8] condensate, and they can constitute all or part of cosmological dark matter [9]. The role of 
i-C supersymmetric Q-balls in cosmology and astrophysics has been the subject of many studies [10-38]. However, it was 
^p-i always assumed that the scalar field that makes up the Q-ball has its VEV aligned exactly along the flat direction. 
O I In this paper we will show that, for a class of flat directions, this is not the case, and we will describe the true ground 
D [ states of the corresponding Q-balls. 

^ • What makes the Q-ball stable is a combination of the baryon number conservation and the energy conservation. 
By construction, a Q-ball has a lower mass than a collection of the scalar particles carrying the same baryon number. 
However, one should also consider the possibility of a decay into fermions [39]. Very large Q-balls are likely to form 
I> from the fragmentation of the AD condensate [9, 36, 40, 41]. For such large Q-balls with the baryon number Qb, the 
)^ ■ mass per unit baryon number decreases with Qs as Q'^^^ . When the mass per unit baryon number is below ra^, the 
proton mass, the Q-ball has a lower mass than a collection of protons and neutrons with the same baryon number. 
Such a Q-ball cannot decay. The same reasoning would not apply to a purely leptonic Q-ball, because the lightest 
fermions carrying the lepton number, neutrinos, have very small masses, one of which can be arbitrarily close to zero. 

However, in supersymmetric models the scalar fields can carry both the baryon number and the lepton number. 
Many flat directions that admit Q-ball solutions have both the baryon and lepton components related to each other 
\ by the requirement that the so called D-terms in the potential vanish. All the previous analyses of supersymmetric 
Q-balls focused on the baryonic squark component, assuming that the slepton fields just follow the squarks. This 
assumption is reasonable, as long as one considers general features of the stable Q-balls. There is an energetic penalty 
for deviations of the leptonic and baryonic components. However, this energetic penalty may not be sufficient to 
$H prevent the emission of a limited number of neutrinos and electrons from a large baryoleptonic Q-ball. Such emission 
can alter the ground state of the relic Q-balls. In particular, it can give the electric charge to an otherwise neutral 
Q-ball.^ Even one unit of electric charge makes a dramatic difference in the experimental signatures of the relic 
Q-balls [18]. For example, the neutral Q-balls can be detected by Super-Kamiokande [18, 42, 43], while this detector 
is not well suited for detection of the charged Q-balls. On the other hand, the charged Q-balls could be detected by 
such detectors as MACRO [18, 44] in some range of parameters. The astrophysical limits on the relic Q-balls may 
depend on the electric charge [17, 27, 28]. It is, therefore, important to understand the ground states of Q-balls for 
the purposes of their detection [18, 42-46]. 

The paper is organized as follows. In section II we review the basic properties of fiat-direction Q-balls in the 
MSSM. In section III we consider the emission of neutrinos from a fiat direction Q-ball associated with the QLd fiat 
direction. In section IV we include the contributions of higher-dimensional operators. In section V we calculate the 
emission of electrons that is energetically allowed. In section VI we show that the relative rates of decay for neutrino 
and electron emission can favor a positively charged Q-ball. 
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^ Electric neutrality is a necessary condition for the cancellation of the D-terms. 
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II. FLAT-DIRECTION Q-BALLS 



Minimal Supersymmetric Standard Model (MSSM) has a large number of flat directions [47]. Let us begin by 
considering a toy model that will capture the main features of Q-balls in the MSSM. 

A flat direction in general can be parameterized by a single scalar degree of freedom, although for many flat directions 
the relation between the gauge eigenstates and the flat direction parameter can be non-trivial [48]. A Q-ball can form 
along a flat direction parameterized by a squark field q and a slepton field L. A Q-ball is a nontopological soliton 
which minimizes the energy under the constraint of a fixed global charge [1-3]. Thus we minimize 



\q' + lW + lL' + l\VL\' + Uiq,L) 



E = j d^x 

while keeping the the baryonic and the leptonic charges separately conserved: 

^B = Yij '^^^ (i^oq* - q*doq) 



QL = Y.j d'x {LdoL* - L*doL) . 



(1) 



(2) 



(3) 



One can solve to the above extremization problem by generalizing the method of Lagrange multipliers as used in 
Ref. [49] to the case of multiple conserved quantum numbers Qa- One looks for a minimum of 



a=L,B 



1 

2i 



d^x {adoa* — a*doa) 



(4) 



We use the same approach as in the single field case [49] to derive the time-dependence of the fields which comprise 
the Q-ball. Namely, we exploit the fact that we can rewrite eq. (4) as 



d^x 
1 



VL\^+U{q,L) 



(5) 



where U = U{q,L) — (1/2)w£L^ — {l/2)uj^q^. Since the first two terms are non-negative definite and no other term 
has an explicit time dependence, one can set them to zero by choosing 



L{r,t) = L{r)e"^^\ 
q{r,i) = q{r)e"^^*. 



(6) 



For definiteness, let us now consider the flat direction for which L = q, with the potential given by U{L,q) = 
Uo + X\L^ — q^\^, where Uq is a constant. The value of A has no effect on the Q-ball solution with wl = i^q = along 
the flat direction, which, in the thick- wall regime [49], can be approximated as follows: 



= i Losm{u;r)/{u;r),r < R 
^ ' 1 Li exp{— mir}, r > R 

air) = I 1osm{ur)/{ur),r < R 
\ gi exp{— mgr}, r > R 



(7) 



Here the constants Li,qi and R are chosen so that they minimize while the solutions are continuous at r = ii. 

We have set the radii of the leptonic and the baryonic components eqiial to each other, which results in the lower 
value of £u- For the flat direction Lq = q^. The exponential tails in eq. (7) give a negligible contribution to most of 
the quantities we compute below; in most cases one can assume that the solution vanishes at r = R. We will use this 
approximation in what follows. 

Using the profiles of eq. (7) in the expression for the leptonic and baryonic charges (2), we get, approximately. 



WL — n-\^0\ , 



Qb = 



(8) 

(9) 
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which must be equal by the Lq = qq constraint. 

The mass of the Q-baU and the frequency lo can be found from the minimization procedure described above: 

M{Ql,Qb) = ^U^^\QL + QBf/\ 



(10) 



u^{Ql, Qb) = V2t,uI'\Ql + Qb)-^'^. 



(11) 



One can also estimate the radius R of the flat-direction Q-baU by neglecting the exponential tail in eq. (7) and setting 
L{R) w 0. This yields 



R{Ql,Qb) = 



oj{Ql,Qb) V2 

eqns. (10)-(12) summarize the general features of the Q-ball formed along a flat direction. 



(12) 



III. EMISSION OF NEUTRAL LEPTONS 

As discussed in the introduction, the flat-dircc:tion Q-balls in MSSM arc stabilized by the baryon number and energy 
conservation, but an emission of some small number of leptons is an interesting possibility. Let us now consider a 
flat-direction Q-ball carrying some lepton number and some baryon number, and decaying into an off-flat-direction 
Q-ball plus some number of fcrmionic leptons. In this subsection we consider only the neutrino emission, so that we 
can ignore Coulomb effects. We label all initial Q-ball values (the FD values) with an i and all final state Q-balls 
with an /. 

We assume that the emission of fermions is a small perturbation from the flat direction case, so that the ansatz 
(/)(r) = sintor /{tor) for both squarks and sleptons still holds. Although away from the direction, the value of mass per 
unit charge in the condensate oj is different from its flat-direction value, the equality of squark and slepton oj's still 
holds for energy reasons. This can be seen by setting iVq = iv, ojl = + in eq. (5) and minimizing with respect to 
Sw. Then the energy is 



£ui,Su = {oj + 5lj)Ql + ujQb 
+ j d^x 

We note that the gradient term for each fleld cancels exactly with the quadratic term hidden in U{q, L), so that 



l\Vq\' + l\Vq\'+Uiq,L) 



(13) 



Su,5w = (w + Su)Ql + loQb + \ j d^x\L^ - 
The integral in the above expression is equal 



I{5uj) 



_x_ r/'^ 



2 / sin(a; -|- Suj)r 



r dr. 



{uj + Suj)r 

We use the definition of the U(l) charges in eq. (2) and take Ql ^ Qb the off-flat-direction state to write 



""Jo 



Next, we evaluate I{6oj) using sin(w -|- Su>)r « {5u)r) cos wr -|- sinwr. Keeping terms to 0{Su)) gives 

J.. , ^XQISCO^ .2 2 , 

lyouj) = g / sm wrcos uirdr 

Jo 

Minimization of fa;,5ai with respect to 6u) yields 



6u = — 



(14) 



(15) 



(16) 



(17) 



(18) 
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For realistically large values oi Qb [9], Soj is very small compared to ui which means Wl ~ Ws- We will neglect the 
difference between w/, and ujb in the remainder of the paper. 

Let us now evaluate the final state (off-flat-direction) Q-ball mass 

Mf = LOfiQf + Qb) + ^ +xj <fx\L^ - gT (19) 

Although this is a (i? + L)-ball, we are interested only in allowing the sleptonic part to change its lepton number. 
Thus to ensure that Qb does not change as we vary L, we demand that 

2 ^2 

% - %■ (20) 

We recall that on the flat direction Li = qi, but this does not hold away from the flat direction. To calculate the 
potential away from the flat direction, we use 

\L'-<ir={^y\L}-,}\'. (21) 
Using the deflnition of the lepton charge, eq. (2), one can relate final and initial leptonic amplitudes: 



Q.) v-.;"* = laJ '''' 

Using the flat-direction expression for slepton amplitude (2) and rewriting in terms of N = Qi — Qf, we obtain a 
simple approximate expression for the potential term: 

xJd'x\L'-q'\' = ^^N\ (23) 

where a4 0.67 is a numerical constant defined by 

r ('J^Yx^dx. (24) 



In particular, a4 — Si{2i:) — \Si{AiT) « 0.67 and « 0.39. Using this expression, we can carry out the minimization 
of eq. (19) and find that the Lagrange multiplier oj is 

^ / 47r^t/o V" ^ ( 1 v'u 1 yi' 

'^^ W{Qf + QB) + Xaim) '^'\l-N/Qi) \l + Xa^N^ /ir^Qf ) ^ ' 

Thus in the limit Qi,Qf ^ N, we see that w/ « uji. We therefore take u)f = in the following. We will verify 
the self-consistency of this assumption ex post facto. Using the value of w/ away from the fiat direction, one can 
compute the final state Q-ball mass as well as the number of neutrinos emitted in the transition from flat-direction to 
ofF-flat-direction states. Since the final state Q-ball is lighter than the original fiat direction Q-ball the excess energy 
is released in neutrinos: 

AM = Mi-Mf= UiN - ^^7V2 (26) 

For the emission to be energetically allowed, one must have AM > Nmi,. This gives an upper bound on the number 
of neutrinos that can be emitted as the Q-ball transitions from the initial flat direction state to a state with the VEV 
that is slightly off the flat direction: 



For A ~ 0.1 and an initial lepton number Qi = 10^"* and u)i ~ Uq^^Q~^^^ ~ 1 MeV, as expected from cosmology [7- 
9], this gives the number of emitted neutrinos A'^ < 460. 

We note in passing that tUv must be less than the uj{Q) for emission to be possible. The physical reason for it is 
clear, since u) can be thought of as the mass of a scalar particle in the condensate. In most models u){Q) » rrii, [7, 8]. 

The potential we have chosen in eq. (19) is of the form one often encounters for MSSM flat directions in gauge 
mediated models [5]. In particular, the results of this section apply directly to the flat directions QdL and QQQL [47]. 
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IV. EXAMPLE: Qid2Li DIRECTION IN THE MSSM 

Let us now analyze a realistic model in which we consider Q-ball formation along the Q\d2Li direction. This 

direction has been studied in detail by Dine, Randall and Thomas [50, 51] in the context of a gravity-mediated Affleck- 
Dine baryogenesis. This direction is useful for our purposes since it carries nonzero baryon and lepton numbers. It is 
useful in the AD scenario because it carries nonzero B — L (in particular B — L = -1), which it must if one would like 
to have a nonzero baryon number after B + L violating sphaleron processes go out of equilibrium. This flat direction 
may be parameterized by the complex field (j) which is the VEV given to the squarks and slepton fields 

where color indices are suppressed and subscripts label generation. The scalar potential on the flat direction is iden- 
tically zero in the supersymmetric limit, since the i^-terms and D-tevms along this direction vanish. This degeneracy 
is lifted however, by soft SUSY breaking terms, as well as by the higher-dimensional operators. 
The SU(2)xU(l) D-terms^ are 

Ud = ^{\Qif~\Lify + ^{\Qi\'-S\L,f + 2\d2f)\ (29) 

where g = e/ sinOw is the SU(2) coupling and g' = e/ cos^w is the U(1)y coupling. When Q\ = L\ = d2 =^ (t> the 
D-terms in the potential vanish as they must. 

The non-renormalizable superpotential along a flat direction has the form 



Wnr = -T^X'' = , 0", (30) 



where n = mk > 4 and X ~ 0™ is a gauge invariant composite operator of the fields that make up the flat direction. 
For the case at hand X = Qid2Li = 0'^. M is a large mass scale usually taken to be either the GUT or the Planck 
scale. This gives a contribution to the scalar potential 

V^nic^) = ^^i<^*^r-'- (31) 
For the flat direction of interest, the lowest order contribution is for n = 4, so 

VNR{ct>)=^-^\<l>\'. (32) 

The additional contribution to AM can affect the decay of the flat direction Q-ball. The new contribution is 
AM D / d'^x{V]srR{(l)FD) — Vnr{4>g)), where 0g and (ppD are the field values in the ground state and flat direction 
configurations, respectively. Since both squark and slepton fleld amplitudes in this case are different from their 
flat-direction values, they both make a contribution to the non-renormalizable part of the potential: 

(AM);vR =%jd'x {\L,f - \Lf\^) + {\q^f - \qf\^) . (33) 
Using the values for Lf and Qf found from the previous section, and assuming N <^Qi, we obtain 

(AM);vH- 247r6M2 " ^^^^ 
This translates into the following bound on N 



Taking A ^ 0.1, M ^ 10^^ GeV, we see that, with the inclusion of the non-renormalizable term, the neutrino 
emission becomes energetically forbidden for Qi > 10^^. 



^ The SU(3) D-terms vanish because the two squark fields, Q and d, have the same amphtudes and the same value of uj. 
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V. CHARGED LEPTON EMISSION FOR QDL 

In addition to the emission of neutrinos, there are, of course, other channels of decay that can alter the Q-ball's 

Icpton number, while leaving the baryon number unchanged. The transition from the flat direction Q-ball to the true 
ground state could be accomplished via the charged lepton emission, for example, via emission of electrons, muons, 
and tau leptons. As before, the emission is only possible as long as lu is greater than the mass of the lepton. In most 
cases, the relic Q-balls have the value of lu below the muon mass [7, 8]. However, the emission of electrons remains 
a possibility. One expects the number of charged leptons emitted from the Q-ball to be smaller than the number of 
emitted neutral leptons because the emission of a each electron adds electrostatic energy to the final state Q-ball. 
The primary difference in this case is that the final state Q-ball is electrically charged. Now the Q-ball mass with 
electric charge eNe is given by 

Mf = u;f{Qf + Q^) + ±^ + + A J d'x\L^ - g^- (36) 

Here we have used the expression for the electrostatic energy of a uniformly charged sphere and R{Q) = n/Lv{Q). As 
before we can find an upper bound on the number of emitted electrons from AM > NrUe- This yields 

N<(^^](l-^). (37) 

Taking A = 0.1 and Ui » rUe, we obtain the bound A'' < 280 on the number of electrons emitted when we consider 
only the electrically charged decay channel 

As before, one must have Lo{Qi) > rrie for emission to be kinematically allowed. Since J/g^^Q^^^^' implies 
Qi < 10^^. If this bound is not satisfied than electron emission is energetically forbidden. We note that the same 
constraint for the neutrino case would allow much larger values: Q, < 10^^ (here we have taken rrii, 0.1 eV). 



VI. COMPETING DECAY CHANNELS 



The relic Q-balls can emit both neutrinos and electrons until they reach the true ground state describe above. In 

each particular case, the history of the AD condensate fragmentation and the subsequent collisions and evaporations 
of Q-balls should be considered. However, it is instructive to consider a simplified history of Q-ball formation. Let us 
assume that at some initial time the Q-ball is in its flat-direction state. This state would be the ground state of an 
infinitely large Q-ball because the mass per unit global charge becomes smaller than the mass of the lightest massive 
fermion in the limit Q ^ oo. The AD condensate, from which the Q-balls form by fragmentation, can be thought 
of as a very large, superhorizon-size Q-ball. Hence, it is reasonable to take the flat direction as the initial state of a 
Q-ball and to consider relaxation of this state into the ground state by emission of fermions, such as neutrinos and 
electrons. The relative rates of decay determine the final electric charge of the Q-ball. Let us find the lowest energy 
state achieved when both decay channels are open. The Coulomb forces disfavor the electron emission. With both 
channels open, the condition for decay AM > J^i NiUii becomes 



207r2 



> Ne,me + N„m^. 



(38) 



With two lepton decay modes, the energetic constraint does not give us a fixed upper bound on N^, but instead 
merely carves out an energetically allowed region in the A^e — N^, plane (See Fig. 1). This region in interesting in that 
it allows the electron emission of a few hundred. Of course, such a larger emission is not realistic since the rate of 
electron emission is expected to be much smaller than the rate of neutrino emission due to the Coulomb suppression. 
It is, therefore, the emission rates that determine the final ground state. 

For a sufficiently large Q-ball the rate of decay into fermions is limited by the number of fermionic states that can 
form inside and cross the boundary per unit time [39] . The maximal rate of fermion emission from the surface of a 
Q-ball [39] is 



dQ 
dt 



' 1927r2 ' 



(39) 
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where A is the Q-baU surface area. This rate is found by calculating the expectation value of the normal component 
(n • j) of the lepton current j in the frequency range < A; < w, as in Ref. [39]'^: 



1 



(27r)3 



.1 .27r 

k dk / cos 6d cos 9 / i 

^0 "'0 



(40) 



This calculation only applies to electrically neutral, massless neutrinos. To extend this result to massive, electrically 
charged leptons (electrons) one must take into account the Coulomb interactions that change the upper bound on the 
frequency. Let us denote this upper bound kmax- For massive, neutral particles k^ax = \/ lo^ — . Including the 
effect of Q-ball charging, we obtain 



2 

max 



207r2 ' 



(41) 



where Zq is the electric charge of the Q-ball (this is N in previous notation). For any type of lepton emission the 
rate of fermion emission from the surface is 



dt 



< 



247^2 '^max- 



(42) 



We can use eq. (42) to track the electron and neutrino emission, in order to determine where in the — plane 
a Q-ball will end. We plot this trajectory in the Ne — N^, plane as the solid curve in Fig. 1. The intersection of the 
energetically allowed boundary (dashed curve) and the dynamical path along which the Q-ball evolves (solid curve) 
gives the maximal electric charge the Q-ball can acquire starting from a (electrically neutral) flat-direction Q-ball. 
This upper bound* is Ne < 25. 



FIG. 1: The region below the dashed curve contains the allowed parameters for the decays from FD Q-balls. The solid curve 
shows the trajectory in the N^-Ne plane determined by the respective decay rates into the electrons and neutrinos. Here we 
show a representative initial Q-ball charge, Q = 10^*. 



^ Our upper limit is lo/2, rather than oj, since the model of Ref. [39] has a scalar decaying into two fermions. 
In the massless limit we find only a slightly different value for the number of emitted electrons A^e < — aV 20/3 : 
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VII. CONCLUSION 

A large number of flat directions in MSSM and other supersymmetric models have both baryon number and lepton 
number simultaneously. Q-balls with VEVs along these flat directions can be entirely stable: they owe their stability 
to the fact that a collection of protons and neutrons with the same baryon number would have a larger mass than 
the Q-ball mass. However, the true ground state of such Q-balls can be different from the naive ground state used in 
the literature. In general, the lepton component of a mixed baryoleptonic Q-ball deviates from the flat direction. We 
have shown this by starting with the naive, flat-direction ground state and following its decay into the true ground 
state plus some number of leptons. The true ground state can be electrically charged, which can have important 
ramifications for the experimental search for relic Q-balls. 

Acknowledgments 

We thank John M. Cornwall for very helpful discussions. We also thank Antoine Calvez and Anupam Mazumdar 
for helpful comments. A.K. thanks Aspen Center for Physics for hospilality. This work was supported in part by the 
DOE grant DE-FG03-91ER40662 and by the NASA ATFP grant NNX08AL48G. 



[1] R. Friedberg, T. D. Lee and A. Sirlin, Pliys. Rev. D13, 2739 (1976). 

[2] S. R. Coleman, Nucl. Phys. B262, 263 (1985). 

[3] T. D. Lcc and Y. Pang, Phys. Rept. 221, 251 (1992). 

[4] A. Kusonko, Phys. Lett. B405, 108 (1997), [hep-ph/9704273]. 

[5] G. R. Dvali, A. Kusenko and M. E. Shaposhnikov, Phys. Lett. B417, 99 (1998), [hep-ph/9707423] . 

[6] L Affleck and M. Dine, Nucl. Phys. B249, 361 (1985). 

[7] K. Enqvist and A. Mazumdar, Pliys. Rept. 380, 99 (2003), [hep-ph/0209244]. 

[8] M. Dine and A. Kusenko, Rev. Mod. Pfiys. 76, 1 (2004), [hcp-ph/0303065]. 

[9] A. Kusonko and M. E. Sliaposlmikov, Phys. Lett. B418, 46 (1998), [hep-ph/9709492]. 

[10] J. A. Fricman, G. B. Gelmini, M. Glciser and E. W. Kolb, Phys. Rev. Lett. 60, 2101 (1988). 

[11] J. A. Fricman, A. V. Olinto, M. Gleiser and C. Alcock, Phys. Rev. D40, 3241 (1989). 

[12] K. Griest and E. W. Kolb, Phys. Rev. D40, 3231 (1989). 

[13] K. Griest, E. W. Kolb and A. Massarotti, Phys. Rev. D40, 3529 (1989). 

[14] K. Enqvist and J. McDonald, Phys. Lett. B425, 309 (1998), [hcp-ph/9711514]. 

[15] K. Enqvist and J. McDonald, Phys. Rev. Lett. 81, 3071 (1998), [hep-ph/9806213]. 

[16] K. Enqvist and J. McDonald, Phys. Lett. B440, 59 (1998), [hep-ph/9807269]. 

[17] A. Kusenko, M. E. Shaposhnikov, P. G. Tinyakov and L I. Tkachev, Phys. Lett. B423, 104 (1998), [hep-ph/9801212]. 

[18] A. Kusenko, V. Kuzmin, M. E. Shaposhnikov and P. G. Tinyakov, Phys. Rev. Lett. 80, 3185 (1998), [hep-ph/9712212]. 

[19] A. Kusenko, Phys. Lett. B406, 26 (1997), [hep-ph/9705361] . 

[20] M. Laine and M. E. Shaposhnikov, Nucl. Phys. B532, 376 (1998), [hcp-ph/9804237]. 

[21] K. Enqvist and J. McDonald, Nucl. Phys. B538, 321 (1999), [hGp-ph/9803380]. 

[22] K. Enqvist and J. McDonald, Phys. Rev. Lett. 83, 2510 (1999), [hcp-ph/9811412]. 

[23] M. Axenides, S. Komineas, L. Perivolaropoulos and M. Floratos, Phys. Rev. D61, 085006 (2000), [hep-ph/9910388]. 

[24] R. Banerjee and K. Jedamzik, Phys. Lett. B484, 278 (2000), [hep-ph/0005031]. 

[25] R. Battye and P. Sutcliffe, Nucl. Phys. B590, 329 (2000), [hep-th/0003252]. 

[26] R. AUahvordi, A. Mazumdar and A. Ozpinoci, Phys. Rev. D65, 125003 (2002), [hcp-ph/0203062]. 

[27] A. Kusenko, L. Lovoridgo and M. Shaposhnikov, Phys. Rev. D72, 025015 (2005), [hcp-ph/0405044]. 

[28] A. Kusenko, L. C. Lovcridgc and M. Shaposhnikov, ,JCAP 0508, Oil (2005), [astro-ph/0507225]. 

[29] M. Berkooz, D. J. H. Chung and T. Volansky, Phys. Rev. Lett. 96, 031303 (2006), [hep-ph/0510186]. 

[30] M. Berkooz, D. J. H. Chung and T. Volansky, Phys. Rev. D73, 063526 (2006), [hep-ph/0507218]. 

[31] A. Kusenko and A. Mazumdar, 0807.4554. 

[32] M. C. Johnson and M. Kamionkowski, 0805.1748. 

[33] S. Kasuya, M. Kawasaki and F. Takahashi, 0805.4245. 

[34] N. Sakai and M. Sasaki, Prog. Tlieor. Phys. 119, 929 (2008), [0712.1450]. 

[35] L. CampaneUi and M. Ruggieri, Phys. Rev. D77, 043504 (2008), [0712.3669]. 

[36] S. Kasuya and M. Kawasaki, Phys. Rev. D62, 023512 (2000), [hep-ph/0002285]. 

[37] M. Kawasaki, K. Kenya and F. Takahashi, Phys. Lett. B619, 233 (2005), [hep-ph/0504105]. 

[38] S. Kasuya and F. Takahashi, JCAP 0711, 019 (2007), [0709.2634]. 

[39] A. G. Cohen, S. R. Coleman, H. Georgi and A. Manohar, Nucl. Phys. B272, 301 (1986). 

[40] S. Kasuya and M. Kawasaki, Phys. Rev. D61, 041301 (2000), [hep-ph/9909509]. 



9 



[41] S. Kasuya and M. Kawasaki, Phys. Rev. Lett. 85, 2677 (2000), [hcp-ph/0006128]. 

[42] J. Arafune, T. Yoshida, S. Nakamura and K. Ogurc, Phys. Rev. D62, 105013 (2000), [hep-ph/0005103]. 
[43] Supor-Kamiokande, Y. Takenaga et al, Phys. Lett. B647, 18 (2007), [hep-ox/0608057]. 
[44] MACRO, M. Ambrosio et al., Eur. Phys. J. C13, 453 (2000), [hep-ex/9904031]. 
[45] L A. Belolaptikov et al., astro-ph/9802223. 
[46] SLIM, S. Cccchini et al, (2008), [0805.1797]. 

[47] T. Gherghctta, C. F. Kolda and S. P. Martin, Nud. Phys. B468, 37 (1996), [hep-ph/9510370]. 
[48] K. Enqvist, A. Jokinen and A. Mazumdar, JCAP 0401, 008 (2004), [hep-ph/0311336]. 
[49] A. Kusenko, Phys. Lett. B404, 285 (1997), [hep-th/9704073]. 

[50] M. Dine, L. Randall and S. D. Thomas, Phys. Rev. Lett. 75, 398 (1995), [hep-ph/9503303]. 
[51] M. Dine, L. Randall and S. D. Thomas, Nucl. Phys. B458, 291 (1996), [hep-ph/9507453] . 



